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Abstract 

The metric on the moduli space of SU(2) charge four BPS monopoles with tetra¬ 
hedral symmetry is calculated using numerical methods. In the asymptotic region, 
in which the four monopoles are located on the vertices of a large tetrahedron, the 
metric is in excellent agreement with the point particle metric. We find that the four 
monopoles are accelerated through the cubic monopole configuration and compute 
the time advance. Numerical evidence is presented for a remarkable equivalence be¬ 
tween a proper distance in the 4-monopole moduli space and a related proper distance 
in the point particle moduli space. This equivalence implies that the approximation 
to the time advance (and WKB quantum phase shift) calculated using the point 
particle derived metric is exact. 
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1 Introduction 


The dynamics of SU(2) BPS monopoles may be approximated by the time evolution of 
a finite number of collective coordinates 1 , 0 . In this moduli space approach the dynamics 
of k monopoles is approximated by geodesic motion on the A;-monopole moduli space Ad*,, 
which is a 4A'-dimensional manifold. To study monopole dynamics therefore requires the 
construction of the metric on AR, which is determined by the kinetic part of the field 
theory action. In the case k = 2, Atiyah and Hitchin |3j were able to calculate the metric 
using indirect methods and making use of its hyperkahler property. However, for k > 2 
the problem is a more difficult one and no metrics have yet been calculated. Recently, 
Gibbons and Manton Jfj have calculated, for general k , the asymptotic metric on regions 
of M.k which describe well-separated monopoles. This asymptotic metric is obtained by 
treating the monopoles as point particles and is of a generalized Taub-NUT form. 

The moduli space A/”, of tetrahedrally symmetric 4-monopoles, is a one-dimensional 
totally geodesic submanifold of Ad 4, and the associated four-monopole scattering process 
has been investigated in some detail [0. The fact that J\f is one-dimensional allows the 
monopole trajectories to be determined even though the metric is not known. In this paper 
we construct the metric on A f by working with Nahm data and using numerical methods. 
In the asymptotic region, in which the four monopoles are located on the vertices of a 
large tetrahedron, the metric is in excellent agreement with the point particle metric. We 
use this metric to calculate the time advance/delay. We also provide numerical evidence 
for the following rather remarkable equivalence. Let l be a good global coordinate on 
A/", such that l = 0 represents coincident monopoles and large l represents well-separated 
monopoles. Then the proper distance from the point l = 0 to a point with / large is equal to 
the proper distance from the singularity to the same point l in the generalized Taub-NUT 
space. This equivalence implies that the approximation to the time advance (and quantum 
phase shift) calculated using the point particle derived metric is exact. This is similar to 
a previous numerical result found in the 2-monopole case 0, and suggests that it may be 
a general feature of the point particle metric. 

2 Four monopoles with tetrahedral symmetry 

In this section we recall the results || on tetrahedrally symmetric charge four monopoles 
that we shall require later. Monopoles are equivalent to various other kinds of mathemat¬ 
ical creatures, and here we shall use two of these; namely, spectral curves and Nahm data. 
Spectral curves 0 are algebraic curves in the holomorphic tangent bundle to the Riemann 
sphere. Let ( be the standard inhomogeneous coordinate on the base space and y the fibre 
coordinate. Then a 4-monopole with tetrahedral symmetry has a spectral curve 

if + i36aK 3 ??C(C 4 - 1) + 3 k 4 (C 8 + 14( 4 + 1) = 0. (2.1) 

where a <G (— a c , a c ), with a c = 3~ 5 / 4 a/ 2, and k is half the real period of the elliptic curve 

y 2 = A(x 3 -x + 3a 2 ). (2.2) 
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Hence there is a one-parameter family of tetrahedrally symmetric 4-monopoles. Since 
this family of monopoles is singled out from the general 4-monopole configuration by the 
imposition of a symmetry, this implies that the corresponding submanifold J\f C Ad 4 
is totally geodesic. The associated monopole dynamics has been studied in detail and 
describes the scattering of four monopoles which arc initially well-separated and positioned 
on the vertices of a contracting regular tetrahedron. As the monopoles merge they scatter 
instantaneously through a configuration with cubic symmetry and emerge on the vertices of 
an expanding tetrahedron dual to the incoming one. For the purposes of constructing the 
metric on J\f we require a good global coordinate which we can identify with the distance 
of each monopole from the origin, at least when the monopoles are well-separated. 

Let x 4 , x 2 , x 3 , x 4 be the four points, each a distance |/| from the origin, given by 


xx = (-/,-/, -0^ 

x 2 = (— 1 , + 1 , +O775 

x 3 = (+/, T/, —(2.3) 
x 4 = (+Z, —l, +0t§- 

They are the vertices of the tetrahedron on which the monopoles are located when they 
are well-separated. For well-separated monopoles the asymptotic spectral curve can be 
obtained as a product of the individual monopole’s spectral curves. If |/| is large and we 
take the four monopolcs to have positions xx,X 2 ,X 3 ,x 4 , then we obtain the asymptotic 
spectral curve 

1 6 4 

^ + ^V(C 4 - 1) + /(C 8 + 14C 4 + 1 ) = o. (2.4) 


The spectral curve (|2.1|) has this form in the limit a —*• a c , upon which k —> oo. By 
comparing ([2 . 1|) and (|2.4|) we see that we can make the identification 


l = A a 1/3 n, where A = 3 7 / e 2~ 2/3 . (2.5) 


At a = 0 (|2T| ) is the spectral curve of the cubic 4-monopolc p, |JJ, which has all four Higgs 
zeros at the origin. If we define, in the usual way, the positions of the monopoles to be 
given by the zeros of the Higgs held, then l = 0 is when all four monopoles have zero 
distance from the origin. Hence l G 1R is a good global coordinate on A f with a natural 
interpretation as the distance of each monopole from the origin. 

We have used the spectral curve approach to monopoles to identify a convenient co¬ 
ordinate on AT, but in order to discuss the metric we now need to turn to the ADHMN 
formulation 0 0. This is usually presented in terms of Nahm data consisting of three 
Nahm matrices (Ti, T 2 , T 3 ), but in order to discuss the metric we must, following Donald¬ 
son 0 , introduce a fourth Nahm matrix Tq. Then we have that charge k monopoles are 


equivalent to Nahm data (To, Tx, T 2 , T 3 ), which are four k x k matrices which depend on a 
real parameter s € [0, 2] and satisfy the following; 

(i) Nahm’s equation 


7/ rj~t -j 

+ [To,n = -£ijfc[Tp TV] 


i = 1,2,3 


( 2 . 6 ) 
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(ii) T 0 is regular for s G [0, 2], Tj(s), i — 1,2, 3, is regular for s G (0, 2) and has simple 
poles at s = 0 and s — 2, 


(iii) the matrix residues of (Tj, T 2 , Tj) at each pole form the irreducible /c-dimensional 
representation of SU(2), 

(iv) 7i(s) = -T) f (s), i = 0,1,2,3, 

(v) T)(s) = T/(2 — s), i = 0,1,2,3. 


Let G be the group of analytic su(/c)-valued functions h(s), for s G [0, 2], which are the 
identity at s = 0 and s = 2, and satisfy /i*(2 — s) = /r _1 (s). Then gauge transformations 
/iGG act on Nahm data as 


T 0 -> hT 0 h~ l 

dh , 

_ ~r h 

ds 

(2.7) 

i ^ hTih - 1 

i = 1,2,3. 

(2.8) 


Note that the gauge To = 0 may always be chosen, which is why this fourth Nahm 
matrix is usually not introduced. However, when discussing the metric on Nahm data we 
need to consider the action of the gauge group and so this extra Nahm matrix needs to be 
kept, at least temporarily. 

In the gauge T 0 = 0 the Nahm data corresponding to a tetrahedrally symmetric 4- 
monopole, whose spectral curve we have discussed above, is given by |JJ 

Ti(s) = x(s)Xi + y(s)Yi +z(s)Zi * = 1,2,3 (2.9) 

where x, y, z are the real functions 

x(s) = 

y(s) = 

z{s) = 

Here p is the Weierstrass function satisfying 

p' 2 = 4p 3 - 4p + 12a 2 (2.13) 


K ( I—. - 7 1 p’(ks) 

- -2 Jp(Ks) + ~ ~ T ~r 
5 V v v 4 p(ks) 


K 
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an 


2p(n s) 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


with prime denoting differentiation with respect to the argument. The tetrahedrally sym¬ 
metric Nahm triplets are 
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Although this Nahm data does not satisfy property (v) this can be achieved by a suitable 
change of basis. 

In the next section we shall use this Nahm data to calculate the metric on J\f. 


3 Calculation of the metric 


It is known that the transformation between the monopole moduli space metric and 

We can therefore calculate the metric, 


the metric on Nahm data is an isometry [12], [13 


g(l), on J\f by computing the metric on the Nahm data given in the previous section. This 
requires the computation of the tangent vector (Vo, Vi, V 2 , V 3 ) corresponding to the point 
with Nahm data (To, Tf, T 2 , T 3 ). In principal, since the Nahm data is explicitly known, this 
could be achieved by direct differentiation, 


Vi = 


dTj 
dl ' 


(3.1) 


However, this is not a practical way to proceed since the Weierstrass function ( |2.13| ), in 
terms of which the Nahm data is given, itself depends on l (through its dependence on a). 
Instead we calculate the tangent space to A f by solving the linearized Nahm equation 


and 


Vi + [Vq, Ti\ + [Tq, 14] — e ijk[Tj, V k \ * — 1,2,3 


Vo + Vi] — 0 

i=0 


(3.2) 


(3.3) 


where Vi, i = 0,1, 2, 3, is an analytic sw(4)-valued function of s G [0,2], Dot denotes 
differentiation with respect to s. The metric on Nahm data is then given by 


g(l) = ~(l j 2 Vi(V?)ds 

J « „0 


(3.4) 


where tr denotes trace and is a normalization constant. 

From now on we use the gauge freedom to set T 0 = 0. Equation (|3.3p is the background 
gauge constraint which ensures that the tangent vectors we compute are horizontal ie 
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that the tangent vectors are orthogonal to the gauge orbits. The tangent vectors are 
tetrahedrally symmetric so we may write 


Vi — Q\X l + q 2 Yi + q 3 Zi i — 1,2,3 


(3.5) 


where q = (qi, q 2 , (fe)* is an analytic real 3-vector function of s G [0, 2], It is easily checked 
that the tetrahedral symmetry of the Nahm triplets implies the following identities 

EPQ, I'.l = EPQ, z,} = E1U, Z,} = 0. (3.6) 

2=1 2=1 2=1 


Substituting ( |2.9| ) and (|3.5| ) into the background gauge equation 
tities (|3.6|) gives the solution Vo = 0. The remaining equations 
equation for the 3-vector q 


and using 
become the 


the iden- 
following 


q = A/q where M 


4x —96 y —12z/5 

—6 y —16?/ — 6a; —6z/5 

—4z —32 z —4x — 32 y 


(3.7) 


Substituting ( |i.5|) into (|3.4| ) gives the metric in terms of q as 


r 2 

g = 120 / (5 ql + 80^ + 3 ql) ds. (3.8) 

Jo 

The ordinary differential equation (|3.7|) has regular-singular points at s = 0 and s = 2, 
since the functions appearing in M have first order poles at these points. Analysis of the 
initial value problem at s = 0 reveals that there is a two-dimensional family of solutions to 
( 0 ) which are normalizable for s G [0,2). They are given by the two-parameter, ai,a 2 , 
family of initial conditions 


q ~ (0, cxis 3 , a 2 s 2 ) t as s ~ 0. (3.9) 

Repeating the analysis for the initial value problem at s = 2 gives a two-parameter, /3±, f3 2 , 
family of normalizable solutions for s € (0, 2], with initial conditions 

q ~ (16/?i(2 — s) 3 , 3/?i(2 — s) 3 , /3 2 (2 — s) 2 ) t as s ~ 2. (3.10) 


The solution we require is the one-parameter family which is normalizable in the closed 
interval s G [0, 2], For later convenience we take a 2 to be the free parameter which describes 
this family of solutions. To compute these solutions by solving an initial value problem 
would be a difficult shooting problem if it were not for the fact that equation (|3.7| ) is linear. 
This reduces the task to a simple problem in linear algebra which we implement as follows. 
Given a value for ck 2 , say a, let Pi(s) denote the solution q(s) of ( |3.7|) corresponding to 
the initial conditions ( |3.9p with (aq,a( 2 ) = (0,cx). This solution is calculated for s G [0,1]. 
Numerically we compute this solution using a fourth order Runge-Kutta method. Let p 2 (s) 
denote a second solution, but this time with initial conditions (ck i,cx 2 ) = (1,0). Similarly, 
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let P 3 (s) and P 4 (s) be the solutions calculated for s G [1,2] obtained from the initial 
conditions fl3.1C| ) with parameter values (/3i,/3 2 ) = (1,0) and (/3 i,/ 3 2 ) = (0,1) respectively. 
Next form the 3x4 matrix 


U 


Pi(l) Pa(l) P)(l) Pi(l) 


(3.11) 


and find the unique solution of the linear matrix equation 

Uw = 0 


(3.12) 


for w = (1, w 2 , W 3 , u> 4 ) f . Numerically this is performed by row reduction of the matrix U 
followed by back substitution. Then the required solution q(s) is given by 


q(s) 


Pi(s) + w 2 P 2 (s) if 0 < s < 1 
— ^ 3 P 3 (s) — w 4 p 4 (s) if 1 < s < 2 


(3.13) 


To summarize, the above procedure consists in integrating ( |3 . 7| ) twice from each end of the 
interval [0, 2] to the centre and then finding a linear combination of these solutions which 
match at the centre. 

We note that in the special case l — 0 (ie a — 0), which corresponds to the cubic 
monopole, the tangent vector may be calculated explicitly in closed form. In this case the 
third component of q decouples from the other two and we have the solution q = (0, 0, g 3 )* 
with 


93 = 


a 2 

k 2 p(ks) 


(3.14) 


where k and p take their values corresponding to a = 0. 

The next issue we confront is to ensure that the tangent vector we compute is dual to 
the coordinate l. This requires the determination of the correct /-dependent normalization 
factor a 2 . In terms of q the equation (|3T|) becomes (qi, q 2 , qY) = ( dx/dl , dy/dl, dz/dl). We 
calculate the correct normalization factor by considering the third component, q 3 = dz/dl , 
in the limit s —> 0. Substituting the asymptotic behaviour of p(ns) as s ~ 0 into the 
expression (|2.12|) for 2 and comparing with the definition of cx 2 given by ( [hOp we obtain 


a 2 s 2 = —(-a,K 3 s 2 ) 
dl K 2 ’ 


which gives 


ot 2 — 


3 1 2 

2A 3 ' 


(3.15) 


There is now only one constant left to determine, which is the overall metric factor 
This is fixed by the requirement that the metric tends to the sum of the monopole masses 
in the limit of infinite separation ie g(l) —> 16(/r as / —> 00 . 
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We now apply the above numerical scheme to calculate the metric. The integral in 
( [T8D is calculated using a standard composite Simpsons rule. The result is displayed in 
Fig. 1 (solid curve) for 0 < l < 18. We see that the metric is a monotonic increasing 
function of l, which implies that the monopoles speed up as they approach each other and 
are accelerated through the cubic monopole configuration l = 0. 

Recently, Gibbons and Manton ||] have calculated the asymptotic metric, on regions of 
M.k which describe well-separated monopoles, by treating the monopoles as point particles. 
In our case of interest the imposed tetrahedral symmetry implies that all four monopoles 
have the same internal phase. This implies that there are no electric charge differences 
between the monopoles so that we can consider pure monopoles not dyons. For k well- 
separated monopoles with positions x, ; , i — 1,.., k, the point particle lagrangian is || 


£ 


2?r H ^ - 2n J2 

i=l 



(3.16) 


For our case of interest k = 4 and the positions are the vertices of the tetrahedron given 
in equation (|2.3|) . Then the above lagrangian becomes 

£ = \g(l)i 2 where g( 0 = 167r ( 1 - y^)- (3T7) 

The point particle metric g(l) is the approximation to the true metric g(l), and the two 
should agree in the large / limit. In Fig. 1 we plot the point particle metric (dashed curve) 
for comparison with the true metric. We see that indeed the two metrics are in excellent 
agreement in the asymptotic limit. This is a useful check not only on the numerics used 
in this paper but also on the point particle approximation applied in Q. Note that the 
metric g has a singularity (ie a point at which its determinant changes sign) at l = \/6. 

Given monopoles which are initially well-separated with l = L , the time taken, t, for 
the monopoles to scatter and reach this separation again may be computed as 


t = A 



where T is the total kinetic energy in the system and A is the proper distance from L to 
the origin 

A = j L \f9Vjdl. (3.18) 

The time advance St, due to the acceleration of the monopoles, is related to the proper 
distance via 


St 


(4Lv£F-A). 


(3.19) 


Note that the asymptotic behaviour of the metric (which is the point particle metric ( |3.17| )) 
means that the time advance has a logarithmic divergence in the limit L —> oo. 





It is interesting to compare the true time advance with that obtained in the point 
particle approximation. In the case of 2-monopoles the point particle metric is the Taub- 
NUT metric with a negative mass parameter [ 141. This also has a singularity, at a finite 
value ro of the radial distance r between the monopoles. Given a point in M 2 corresponding 
to large r, the Atiyah-Hitchin metric can be used to calculate the proper distance of this 
point to the bolt (which describes coincident monopoles). The proper distance from the 
same point r in Taub-NUT space to the singularity can also be calculated and the two 
compared. It is curious that numerically these two distances are found to agree 0. This 
equivalence, for which there is at present no explanation, implies that the approximation 
to the time advance calculated using the Taub-NUT metric is exact. 

Given that we have the metric g and its point particle approximation g we can investi¬ 
gate the possibility that a similar result to that above exists in the 4-monopole case. The 
proper distance from the singularity in the generalized Taub-NUT space is 


A = [ \fg{V) dl = 4^(7 L(L - V6) + log(^— = )) (3-20) 

V2 7I +v /l-v/6 

~ 4Vtt(L- y|(i og (i|) + i) as L —> oo. 

It is this which we wish to compare with A in the large L limit. Numerically we have 
computed the metric g[l) at values up to l = 18. It can be seen from Fig. 1 that this 
is a reasonably large value since the metrics are very similar for l > 10. Setting L = 18 
in (|3.20|) gives the result A = 89.7. The integration to calculate A at L = 18 from the 
numerical values for y fg is performed with the use of the numerical routines FITPACK. We 
fit a spline under tension to the data values and integrate the resulting spline to obtain the 
result A = 89.2. So, to within the numerical accuracy of the calculation, we find that the 
two answers agree. This implies that there is no relative WKB quantum phase shift from 
the point particle approximation in the quantized dynamics of the above classical motion. 

The result in the two monopole case and the numerical evidence presented here suggests 
that perhaps this feature of the point particle metric is a general one. At present there is 
little explanation for this possible equivalence, but perhaps the answer lies in some global 
geometrical properties of the metrics and the fact that the point particle metric inherits 
the hyperkahler property of the true metric on Aik- 


4 Conclusion 

We have introduced a numerical scheme to calculate the monopole moduli space metric 
from Nahm data. This scheme has been used to calculate the metric on a totally geodesic 
submanifold of the 4-monopole moduli space, corresponding to tetrahedrally symmetric 
monopoles. The results compare well with the asymptotic point particle metric and we 
have presented evidence for a curious exact result using an approximate metric. The scheme 
can be applied to calculate the metric on submanifolds of Mk for which the Nahm data 
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is known. A suitable candidate is the submanifold of Ad 3 obtained by imposing a twisted 
line symmetry on three monopoles [l5j . In this scattering process it appears that the zeros 


of the Higgs field stick at the origin for a finite time interval. A calculation of the metric 
would reveal the time scale over which this sticking takes place. 
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Figure 1: The metric g(l) (solid curve) and the point particle metric g(l) (dashed curve). 
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